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We show that superlattices based on zero-gap semiconductors such as graphene and mercury 
telluride exhibit characteristic Bloch-Zener oscillations that emerge from the coherent superposi- 
tion of Bloch oscillations and multiple Zener tunneling between the electron and hole branch. We 
demonstrate this mechanism by means of wave packet dynamics in various spatially periodically 
modulated nanoribbons subject to an external bias field. The associated Bloch frequencies exhibit 
a peculiar periodic bias dependence which we explain within a two-band model. Supported by 
extensive numerical transport calculations, we show that this effect gives rise to distinct current 
oscillations observable in the I-V characteristics of graphene and mercury telluride superlattices. 

PACS numbers: 72.80.Vp, 73.21. Cd, 85.35.Ds, 85.75.Mm 



I. INTRODUCTION 

Bloch oscillations, the periodic motion of particles in 
a superlattice subject to a constant external field, rep- 
resent a fundamental phenomenon in transport through 
periodic potentials. Predicted already in the early days 
of quantum mechanics^ Bloch oscillations have been 
observed in various fields of physics, ranging from ear- 
lier experiments in semiconductor superlattices^^ via 
cold atoms in optical lattices^'^ to classical optical^'^ 
and acoustic^^ waves. While many aspects of conven- 
tional Bloch oscillations can be explained by a single 
band description, particularly interesting effects arise in 
the case of two coupled minibands^^ energetically sepa- 
rated from further bands. Then partial Zener tunneling 
at avoided crossings of the two minibands can lead to a 
coherent superposition of Bloch oscillations^^'^^, i.e. to 
a splitting, followed by as subsequent recombination of 
a Bloch oscillating wave packet. This gives rise to a va- 
riety of Rabi-type interference phenomena, in particular 
double-periodic motions coined Bloch-Zener (BZ) oscil- 
lations^^^^. Signatures of this effect have already been 
detected in the THz emission of AlGaAs superlatices^^, 
and even the population dynamics have been measured 
recently for light^^ and atomic matter waves^^ in espe- 
cially tailored binary lattices. 

However, materials with a linear Dirac spectrum^^ nat- 
urally serve the effect, since only a small gap is opened by 
a spatially periodic modulation allowing for Zener tun- 
neling between electron and hole states. Such materi- 
als are now at hand with the discovery of graphene^ ^'^^ 
and the advent of topological insulators^^^^ first realized 
in two-dimensional mercury teluride (HgTe) heterostruc- 
tures^^'^^. Interesting phenomena for graphene based 
periodic superstructures have already been theoretically 
predicted like the formation of extra Dirac cones^^^^ and 
the appearance of a negative differential conductance^^. 
Furthermore, recent experiments have realized graphene 
superlattices with periodicities down to a few nm^^. 

This raises the question for the existence of peculiari- 
ties of Bloch oscillations in graphene and topological in- 
sulator superlattices that we address in this manuscript. 



We are not aware of work showing unconventional fea- 
tures in graphene-based Bloch oscillations. Up to now, 
only the semiclassical approach was adapted to a linear 
dispersion^^ and, without reference to Bloch oscillations 
numerical evidence for a negative differential conduc- 
tance was reported^^. We show that besides conventional 
Bloch oscillations, multiple Zener tunneling between the 
coupled electron and hole branches leads to distinct BZ 
oscillations that appear to be naturally present in super- 
lattices made of systems with Dirac-like dispersion. 

This paper is structured as follows: In Section II we 
show the influence of BZ oscillations on the wave packet 
motion in a graphene nanoribbon and the influence on 
the frequency spectrum. Subsequently, we introduce in 
Section III a two-band model to explain the effect in the 
frequency spectrum and the influence of BZ oscillations 
on the electron-hole polarization. In Section IV we show, 
that the occurrence of BZ oscillations can be seen as dis- 
tinct features in the current through graphene nanorib- 
bons. In Section V we present results that feature the 
special frequency pattern of BZ oscillations, as well as 
their signatures in transport, in nanoribbons made of the 
topological insulator mercury telluride. 



II. WAVE-PACKET MOTION IN GRAPHENE 
SUPERLATTICES 

An insight into the dynamics of BZ oscillations can 
be gained by the time-evolution of a wave packet on 
a graphene nanoribbon in presence of a periodic mass 
potential M{x) and a linear electrostatic drift potential 
V{x) as sketched in Fig. 1(a). To this end we model the 
electronic structure of graphene by a conventional tight- 
binding Hamiltonian^^ 

^tb = ^ t cl_^cj,^ + V cl^c,,^ p c{^c,^p (1) 

{i J>,/3 

where {i j) denotes neighbouring unit cells and /3 = ±1 
the sublattice degree of freedom. The initial wave packet 
is created by diagonalizing the periodic Hamiltonian 
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FIG. 1. {Color online) Exemplary setup for Bloch-Zener os- 
cillations in a graphene nanoribbon. a) Sketch of a Gaus- 
sian wave packet in presence of a periodic mass potential 
M{x) = Mo sin(27rx/a) and an electrostatic drift poten- 
tial V{x) = —cEdx. b) Bandstructure of the super lattice 
with small avoided crossing at /c = (nanoribbon width 
W = 10 ao, a = lOVS ao, Mo = 0.1 t). Thick and dashed 
lines show the first and second Bloch band from the metallic 
armchair mode. The gray lines represent higher transversal 
modes. 

H{k) of one unit-cell of the infinite ribbon. By means 
of the transversal eigenfunctions Xniv^k) we create an 
initial electron-like wave packet 

/CO 
X„(y,fc)e"=-e-5'='^'dfc (2) 
-OO 

with a Gaussian broadening S. Since the armchair bound- 
ary mixes the two graphene valleys the wavefunction 
comprises several nodes in lateral direction. The time- 
evolution is calculated by an expansion of the time- 
evolution operator in Chebychev polynomials^^. In pres- 
ence of a periodic mass potential 

M{x) = Mo sin(27rx/a), (3) 

where Mq is the strength of the periodic mass, and a 
the periodicity length, the bandstructure of the superlat- 
tice exhibits a small anti-crossing at /c = and a large 
bandgap between the first and the higher Bloch bands 
as shown in Fig. 1(b). In presence of the a linear drift 
potential 

V{x) = -eEox, (4) 

with Ed as the strength of the drift field, the wave packet 
starts to accelerate. Because of its extent in the longitu- 
dinal direction, the wave packet is localized in momen- 
tum space with a distinct average momentum k{t) in x- 
direction. Given the periodicity of the bandstructure, a 
sawtooth behavior of k{t) is obtained known as Bloch 
oscillations. 

However, the dynamics in a graphene nanoribbon 
shows additional features due to the strong coupling be- 
tween electron and hole states. Therefore we study a 
typical trajectory 

x{t) = {i^{t)\x\m) (5) 

of the center-of-mass (COM) as shown in Fig. 2(a). 
Initially, the wave packet is chosen electron like and 
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FIG. 2. {Color online) Snapshots of a wave packet in the 
course of Bloch-Zener oscillations, a) Center-of-mass mo- 
tion of a wavepacket on a graphene nanoribbon {W — lOao, 
a — 10A/3ao, Mo = O.lt). b) Snapshots of the probability 
distribution of the wavepacket for the corresponding times 
marked with crosses in panel (a). Please note the video of 
the dynamics in the online version of the Supplemental Ma- 
terial^'^. 



a snapshot of the probability distribution is shown in 
Fig. 2(b0). During the first Bloch cycle the probability 
distribution is predominantly to the right of the initial 
position [see Fig. 2(bl)]. This region features a nega- 
tive electrostatic potential and accordingly the part of 
the wave packet with electron character performs Bloch 
oscillations in this region. In Fig. 2(b2) the electron and 
hole parts meet again in momentum space and as a con- 
sequence tunneling from the electron to the hole branch 
is possible as sketched by the bullets in Fig. 1(b). As a re- 
sult the hole-like part of the wave-packet increases and in 
subsequent time steps the COM trajectory reaches neg- 
ative values. The corresponding snapshot at the turning 
point of the the wavefunction in Fig. 2(b3) shows a big 
hole-like state on the left side and a smaller electron- 
like state on the right side. After the next tunneling the 
probability distribution between electron and hole states 
is almost equal, thus the COM motion is strongly sup- 
pressed. Because of the periodic mass potential, the gaps 
between the first Bloch band and higher bands is bigger 
than the gap between electron and hole states as shown in 
Fig. 1(b). As a result the tunneling into higher bands is 
very unlikely and there is no damping of the oscillations 
due to leakage into higher bands. 

To study the dynamics of the tunneling between the 
electron and the hole branch in more detail we perform a 
frequency analysis of the COM motion for different fields 
Ed. The Fourier amplitudes of the dominant frequency 
contributions are visualized by dark colors in Fig. 3. Be- 
sides the conventional Bloch frequency marked by a white 
dashed line, the resulting spectrum shows a pronounced 
interweaving pattern around half of this frequency (black 
dashed line). A stronger periodic potential, and thereby 
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FIG. 3. {Color online) Frequency spectra E = fvuj from the 
center-of-mass motion of a wave packet for varying drift field 
Ed for (a) moderate (Mo = O.lt) and (b) stronger periodic 
potential (Mo = 0.2t). Dark colors represent strong intensi- 
ties. The dashed lines correspond to {1/2, 1, 3/2} times the 
conventional Bloch frequency. 
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FIG. 4. {Color online) a) Bandstructure of the Dirac model 
Hamiltonian (6) for ^; = 1, ^ = 1, a = 1/10, g = 1/2. The 
shaded (yellow) area denotes the integral (9). b) Frequency 
spectrum of the Bloch oscillations for different drift accelera- 
tions a = eEo/h. Solid lines show the frequencies nuj^ +muj~ 
given by Eq. (16), dotted (dashed) lines the strong (weak) 
tunneling limit. 



an increased gap between electron and hole branch, leads 
to a rhombic structure as shown in Fig. 3(b). These pe- 
riodic features in the frequency spectrum arise from the 
interplay between Bloch oscillations and the splitting of 
the wave packet into the electron and hole branches at 
k{t) = (see Fig. 1(a)). The persistent sequence of tun- 
neling events between the two branches and the subse- 
quent interference leads to a new set of frequencies which 
can be unterstood by means of the following model. 



III. ANALYTICAL MODEL FOR BZ 
OSCILLATIONS 

In the following we quantitatively explain these charac- 
teristic BZ features using a periodically modulated one- 
dimensional Dirac model Hamiltonian, 



H{t) = ^ sin {^)a^ 



(6) 



Here a is the period, v is the Fermi velocity and g the 
energy gap between the electron and the hole states. The 
resulting bandstructure is given by 



2{hvla)'^[l-cos{ak)] 



(7) 



as shown in Fig. 4(a). A comparison with the full 
tight-binding calculation of the graphene nanoribbon in 
Fig. 1(b) shows a very good correspondence. An external 
electric drift field Ed enters the equations of motion for 
the quasi-momentum k{t) as hdtk{t) = eEjj leading to a 
time evolution of k{t) = at linear in t where a = eE^/h. 
Conventional Bloch oscillations with frequency ujb = (^a 
arise from the periodicity of k{t) in momentum space 
in the interval [— ^, J]. The phase (j) between the two 
branches accumulated during one oscillation is given by 
a free propagation and thus 



A 

eE 



16v 



D 



(8) 



with 



/Tv/a 
{e^ -e-)dk 
-TV /a 



(9) 



the area in momentum space as depicted in Fig. 4(a). 
This free propagation can be expressed by the matrix 



>i0/2 



e- 



-iS/2 



(10) 



Additional to conventional Bloch oscillations on either 
branch, there is a strong periodic tunneling between the 
electron and the hole states close to the anti-crossing at 
k = 0. There, the Hamiltonian (6) can be linearized 
[dashed lines in Fig. 4(a)], leading to a typical Landau- 
Zener tunneling problem^^^^: 



_ hv at g 
^ g -hv at 



(11) 



The scattering between the different branches is de- 
scribed by 



^0 



(12) 



and 



with the tunneling rate q = 1 — e , S — 2ffiva ' 
^ = J-|-arg(l — i(5)-l-(5(log(5— 1) is an additional tunneling 



phase. From this we can deduce the scattering matrix 



s 



pi(?-0/2) 



(13) 



which describes the time-evolution of the electron and 
hole branch for one Bloch cycle. Using this matrix we 
derive scattering eigenstates 



,± = ^ Vgcos(,^/2 -i)±J\- qsm\<j)/2 - 



0/2 



(14) 
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FIG. 5. {Color online) a) Polarization dependence on the 
phase difference between electron and hole branch of the 
scattering eigenstate (solid line shows the upper spinor 
entry, dashed line the lower spinor entry). Panels (b), 
(c) show the cent er-of- mass motion of an initially electron- 
polarized wavepacket on a graphene nanoribbon superlattice 
M(x) = Mo + V{x) for different drift fields Ed (Mq = 50 meV, 
V{x) = 300meVsin(27rx/a), a = lOVSao). Blue dots de- 
pict regions with a negative amplitude, corresponding to a 
wavepacket with strong hole character. 

with the corresponding eigenvalues e^^^ where 

= arccos ^±^1 - g sin^(^/2 - 0^ • (15) 

The phases f3^ of the scattering eigenstates depend peri- 
odically on the phase difference (j) between electron and 
hole branch. This periodicity leads to two new Bloch 
frequencies 



aa 

TT 



arccos [ ± y/qsm{(j)/2 — ^)] . 



(16) 



Unlike conventional Bloch oscillations these frequencies 
do not simply depend linearly on the drift strength a, 
but show a rapid interweaving pattern strongly changing 
with a, as shown in Fig. 4(b), owing to coherences from 
combined dynamics on the hole and electron branch. The 
limiting cases can be understood as follows. For strong 
coupling, the tunneling rate q ^ leads to a frequency 
cj^ ujb/2 [dotted line in Fig. 4(b)], since for every 
Bloch cycle the states tunnel completely between the two 
branches in momentum space and hence the complete cy- 
cle in position space is twice as long. In the opposite, 
weak coupling limit uj^ ^ aa[l/2± (0/2 — 0/'^] ^ 
leading to a rhombic frequency pattern shown as dashed 
lines in Fig. 4(b). For intermediate tunneling rates the 
frequencies show a smooth transition between these lim- 
iting cases and are in very good agreement with the nu- 
merically calculated spectra of Fig. 3(a,b). 

Furthermore, the scattering eigenstates show a strong 
polarization dependence (electron or hole type character) 
on the phase 0, Eq (8). If the one-dimensional model 
Hamiltonian (6) is considered for g'^ ^ 2h^va^ the tun- 
neling rate g ^ 1, which results in strongly electron or 
hole polarized states foi" almost all values of (j) . The 
absolute value of the spinor entries is always very close to 
one or zero as shown in Fig. 5(a). Nevertheless, the po- 
larization breaks down whenever the difference between 



the phase of the electron and hole branch is 
(^ = 2(n7r + 



- TT 



(17) 



where n G N. This alternating weight of the spinor be- 
tween the electron and hole type states for different drift 
fields Ed can be also deduced from the COM motion of 
wave packets with fixed initial polarization. If the drift 
field is adjusted such that the phase condition (17) is ap- 
proximately satisfied, the COM motion of the initially 
electron-like configuration exhibits oscillations ranging 
from —15 nm to 15 nm for Ejj = 4.61 mV/nm as shown in 
Fig. 5(b). Since conventional Bloch oscillations in a sin- 
gle band are restricted to positive or negative values the 
trajectories imply strong tunneling between the electron 
and hole states. For values of Ed where condition (17) 
is not fulfilled, e.g. Ed — 4.62mV/nm in Fig. 5(c), the 
trajectories of the different polarizations do not signifi- 
cantly cross the origin, thus they preserve their electron- 
hole character. As a consequence, if charge transport 
through a system comprises a transition from electron to 
hole states the current should strongly depend on the BZ 
oscillations within the superstructure. 



IV. TRANSPORT IN GRAPHENE-BASED 
SUPERLATTICES 

In the following, we consider charge transport through 
graphene nanoribbon based superlattices and demon- 
strate that BZ oscillations lead to clear-cut features in 
the I-V characteristics. To this end we model a graphene 
nanoribbon of width W and length L by the tight-binding 
Hamiltoninan of Eq. (1), now with a periodic electro- 
static potential sin(27rx/a) leading to a superlattice 
mini-bandstructure as shown in the inset of Fig. 6(a). 
A small constant mass term M{x) = Mq is additionally 
considered which opens up a gap commonly present in ex- 
periments on graphene nanoribbons^^. We assume a lin- 
ear potential drop eFsD^/^ due to the source-drain volt- 
age FsD between the graphene leads at x = ±L/2. The 
current is calculated by means of the Landauer-Biittiker 
formalism^^ , 



SD 



)=-r 

' h J., 



TiE, Vsd)[/+(^) - f-iE)]dE, (18) 



with the Fermi functions f^{E) = {1 + exp[(i? =F 
ysD/2)/fcBT]}-i. 

As shown in Fig. 6, the current through the nanoribbon 
is governed by a conventional increase with the bias win- 
dow for small FsD, followed by a region of negative dif- 
ferential conductance typical for superlattices. At higher 
bias, FsD > 0.3 V, we observe the emergence of distinct 
current oscillations that get more pronounced with in- 
creasing gap size, see Fig 6(b). Due to the bias be- 
tween source and drain electrode the particles traversing 
the superlattice must change their electron-hole charac- 
ter. However, states performing BZ oscillations exhibit 
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FIG. 6. {Color online) Current- volt age characteristics for 
graphene nanoribbon superlattices (L = 3000A/3ao, W — 
lOao, a — SOv^tto, Vb = 500 meV) for (a) different Fermi 
energies (Mo = 20meV, T = 20K) and (b) different tempera- 
tures (Mo = 50meV, Ep — ^) showing pronounced signatures 
of Bloch-Zener oscillations at higher bias. Arrows mark ex- 
pected peak positions from phase condition (17). Upper inset: 
Bandstructure (for Mo = 20meV), lower inset: Transmission 
map T(^, Vsd) used in Eq. (18) to get the current of panel 
(b); dark colors represent high transmissions. 



transitions between the two carrier types only for cer- 
tain Vsd = eEjjL when the phase (j) fulfills the condition 
of Eq. (17) as shown in the previous section. In conse- 
quence the current is strongly enhanced if this is fulfilled. 
As shown in Fig. 6(a,b) the current peaks calculated by 
Eq. (18) perfectly coincide with the expected voltages 
(marked by vertical arrows) deduced by extracting the 
area A in momentum space from the minibandes around 
the Fermi energy shown as shaded area in the inset of 
Fig. 6(a). Vice versa, the experimental observation of BZ 
peaks in the I-V characteristics would allow for 'measur- 
ing' the miniband structure. 

A closer look at the transmission values T{E^ Vsb) [see 
inset Fig. 6(b)] reveals a rhombic structure which features 
pronounced transmission maxima piled up at these par- 
ticular values of Vsb (dashed lines) . Since these maxima 
are present for various energies in the conductance win- 
dow, the resulting current is fairly independent of the 
exact Fermi energy [see Fig. 6(a)] and temperature [see 
Fig. 6(b)]. 



V. BZ OSCILLATIONS AND TRANSPORT IN 
HGTE-BASED SUPERLATTICES 

A different setup featuring BZ oscillations can be cre- 
ated from a strip etched out of the two-dimensional topo- 
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FIG. 7. (Color online) Bloch and Bloch-Zener oscillations in 
spatially modulated two-dimensional HgTe nanoribbons. a) 
Bandstructure for a HgTe nanoribbon with periodically mod- 
ulated width W{x) (Eq. (20)) ranging from Wo = 300 nm 
to Wi = 50 nm and periodicity a = 200 nm. b) Frequency 
spectrum E = huu of the wave packet center-of-mass motion 
as a function of drift field Ed- Dashed lines indicate the 
frequencies of the Bloch oscillations, c) /-Vsd characteris- 
tics of a nanoribbon with constant width VK = 150 nm and 
electrostatic modulation V{x) = Vb sin(27rx/a). Small verti- 
cal arrows mark expected maxima from phase condition (17). 
Inset: corresponding miniband structure. 



logical insulator based on mercury teluride (HgTe)^^'^^. 
We describe the electronic properties of the underlying 
HgTe heterostructure by the Hamiltonian^^ 



H 



fCk + Mk Ak+ 

Ak_ Ck-Mk 

Ck^Mk -Ak- 

V -Ak+ Ck-Mu 



(19) 



where k± = kx ^ iky, k = k^ -\- ky, Ck = —Dk and 
Mk = M - Bk^. We assume an HgTe/HgCdTe het- 
erostructure with a quantum well width of 7.0 nm fea- 
turing topological edge states and leading to material pa- 
rameters A, M as typically used in literature^^. As 
for the graphene nanoribbon we can create two different 
types of superlattices with a mass-like modulation and 
an electrostatic modulation. For a HgTe strip the mass 
modulation can be achieved by modulating the width of 
the ribbon, for example by 



W{x) = 



Wo + Wi Wo 



-Wi . f27rx\ 



where Wo and Wi are the maximum and minimum width 
and a is the periodicity. The finite width of the HgTe 
nanoribbon allows for a hybridization of the edge states 
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with the same spin at the opposite boundaries leading to 
a smah gap in the bandstructure^^'^^. Accordingly, the 
modulation of the width corresponds to a modulation of 
the mass gap. The resulting miniband structure from the 
two-dimensional system, shown in Fig. 7(a), is obtained 
numerically by Lanczos diagonalization and exhibits var- 
ious Landau- Zener anticrossings within the bulk bandgap 
of HgTe which suggest BZ oscillations. 

In order to study the electron dynamics we calculate 
the COM motion of Gaussian shaped edge-state wave 
packets. Initially, the wave packet is localized on one 
edge and the direction of motion is determined by its 
spin. The array of multiple constrictions enables tunnel- 
ing between the edges. As a consequence an inversion 
of the direction of motion is possible, leading to Bloch 
and BZ oscillations. As shown in Fig. 7(b), the result- 
ing frequency spectrum features the expected rhombic 
pattern in between the frequencies of the conventional 
Bloch oscillations (white dashed lines). Compared to the 
graphene system [see Fig. 3(e)] we observe more com- 
plicated, superimposed structures because of the whole 
sequence of multiple anticrossings in the band structure 
that affect BZ oscillations. 

As for graphene we further study the transport prop- 
erties of HgTe strips of constant width and a periodically 
modulated electrostatic potential resulting in a supercell 
bandstructure shown in the inset of Fig. 7(c). The small 
gap between the electron and the hole states is attributed 
to the finite ribbon width of 150 nm. We chose the Fermi 
energy close to the band crossing of the topological edge 
states and calculate the current using Eq. (18). Besides 
a strong negative-differential conductance at lower bias 
we get the signatures of BZ oscillations for Vsb > 9 mV 
as shown in Fig. 7(c). Similar to the calculations for the 
graphene superlattice the oscillations are independent of 
the exact choice of the Fermi level. The peak positions 
are in good accordance with the expected series of drift 
voltages marked by arrows in Fig. 7(c) obtained from 
Eq. (17), where A is extracted from the bands around 
the Fermi-energy shown as shaded area in the inset. 



VI. CONCLUSION 

In this manuscript we showed that Bloch-Zener oscil- 
lations appear naturally in superlattices made of materi- 
als with a Dirac-like spectrum, highlighting interference 



between electron and hole states. The characteristics 
of these oscillations are explained by a one-dimensional 
model Hamiltonian and numerically confirmed for re- 
alistic setups by means of wave packet simulations for 
graphene and topological insulator ribbons. Further- 
more, we demonstrated that Bloch-Zener oscillations 
manifest themselves as regular sequence of pronounced 
current peaks in quantum transport, besides the well 
know negative differential conductance at low bias, a sig- 
nature of conventional Bloch oscillations. The sequence 
of current peaks associated with the Bloch-Zener oscil- 
lations are intimately linked to the underlying miniband 
structure. 

We suggest transport measurements through graphene 
nanoribbons and HgTe strips as promising experimental 
setups that feature Bloch-Zener oscillations. For single 
layer graphene and topological insulators, the periodic 
electrostatic potential can be imprinted by an array of 
top gates. The gap between the electron and hole states 
can be tuned by the width of the considered nanorib- 
bons. In case of bilayer graphene the gap can also be 
created via a potential difference in z-direction induced 
by top gating. The calculations presented here have 
been performed for clean, disorder free and coherent sys- 
tems. However, preliminary numerical calculations for 
graphene-based superlattices with disorder indicate that 
Bloch-Zener oscillations are still visible if the mean-free- 
path exceeds several periods of the superlattice. This is 
promising with respect to their experimental detection in 
sold-state based samples. 

We finally note that signatures of the Bloch-Zener os- 
cillations presented have been recently observed with ul- 
tracold, fermionic K-atoms due to the Dirac points with 
small mass gaps emerging in tunable optical honeycomb 
lattices^^. 

Note added in proof. Recently, we became aware of Ref. 
47 where the Bloch-Zener oscillations of collective excita- 
tions in narrow zigzag-shaped optical lattices is studied 
theoretically. 
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